Spectroscopy of the unbound: One-neutron knockout reaction of Be 
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Full Faddeev-type calculations are performed for one-neutron knockout reaction of ^*Be on proton 
target at 69 MeV/u incident energy. Inclusive transverse momentum distributions for the outgo- 
ing (^^Be + n) system and semi-inclusive cross sections are presented. A significant proton-core 
single scattering contribution emerges where the valence neutron has non zero angular momentum 
relative to the core. This indicates that distorted-wave impulse approximation is inadequate and 
the complete multiple scattering series must be taken into account for the considered reaction. The 
magnitude of the semi-inclusive cross section at quasifree scattering conditions is a clear signature 
of the angular momentum of the valence nucleon. 
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I. INTRODUCTION 

Knockout-reactions constitute a sensitive tool to inves- 
tigate the single-particle or cluster structure of nuclei. In 
the case of one-neutron knockout, the measurement of 
the momentum distribution of the neutron and the A-1 
system allows to identify, through the shape of the dis- 
tributions, the orbital angular momentum of the struck 
particle. 

Recently one-neutron knockout in the ^''Be scatter- 
ing from a proton target at 69 MeV/u was measured at 
RIKEN in order to obtain information on the unbound 
^^Be [l| and the data analysed through the invariant mass 
method. This is a very powerful tool to study unbound 
nuclear states where the invariant mass is determined by 
the momentum vectors of the outgoing particles the 
advantage of the method is a very good energy resolution. 
Nevertheless, very little is known about the ^^Be system 
and contradictory results have been found both from the 
theoretical and experimental sides. Although it would 
be more appropriate to describe ^'*Be-p scattering in a 
4-body model (^^Be,n,n,p), the exact 4-body scattering 
equations at present can only be solved for the four nu- 
cleon system below breakup threshold. Therefore in this 
paper we use the 3-body Faddeev/Alt, Grassberger, and 
Sandhas (AGS) scattering formalism 0, 0, We ex- 
plore a number of ^^Be configurations and calculate the 
breakup observables to shed light on the structure of the 
unbound system and to stimulate further experimental 
work. 

In the absence of exact many-body reaction calcula- 
tions, Distorted Wave Impulse Approximation (DWIA) 
methods have been applied to study such reactions Q. 
In this approach one assumes that the incoming parti- 
cle collides with the struck nucleon as if it was free and 
the relative motion of the particles in the entrance and 
exit channels is described by distorted waves. Moreover 
further approximations are usually made in standard ap- 
plications of DWIA when evaluating the transition ampli- 



tude for the knockout scattering process A{a, ah)B such 
as: (i) the potential approximation in the entrance chan- 
nel VaA — Vab ~ VaBi'i^aA)', (ii) the factorizatiou approxi- 
mation, which is only exact in Plane Wave Impulse Ap- 
proximation (PWIA); (iii) the on-shell approximation of 
the transition amplitude. Although the validity of some 
of these approximations has not been properly investi- 
gated yet, it has been shown Q that DWIA leads to 
an incomplete and truncated multiple scattering expan- 
sion that is responsible for inaccurate results of ^^Be-p 
breakup observables at intermediate energies. The work 
in Refs. 0, H, Q is based on exact Faddeev/AGS theory 
0, 0, Q and provides not only a benchmark to study ap- 
proximate methods commonly used in nuclear reactions, 
but also a means to obtain a better description of the 
dynamics that drives knockout reactions involving light 
bound systems. 



II. THE FADDEEV/AGS EQUATIONS 

Let us consider 3 different particles denoted 1,2,3, in- 
teracting by means of two-body potentials. The Fad- 
deev/AGS multiple scattering framework is a three-body 
scattering formalism that treats all open channels (elas- 
tic, inelastic, transfer and breakup) on equal footing. 
It is therefore adequate to describe the scattering of 
light bound systems such as halo nuclei where breakup 
thresholds are close to the ground state. In the Fad- 
deev/AGS framework the full scattering amplitude may 
be expressed as a multiple scattering expansion in terms 
of the two-body transition operators for each interacting 
pair. We summarize here the main expressions of the 
formalism, using the odd man out notation appropriate 
for 3-body problems which means, for example, that the 
interaction between the pair (23) is denoted as vi. As- 
suming that the system is non-relativistic, one writes the 
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total Hamiltonian as 

H = Ho + Y,y7^ (1) 

7 

with the kinetic energy operator Hq and the interaction 
for the pair 7. The Hamiltonian can be rewritten as 

H^H^ + V , (2) 

where Ha is the Hamiltonian for channel a 

Ha=Ho+Va, (3) 

and V" represents the sum of interactions external to 
partition a 

^" = E "7 • (4) 

The a = partition corresponds to three free particles 
in the continuum where V'^ is the sum of all pair interac- 
tions. Let us consider the scattering from the initial state 
a to the final state /3. The operators U^°', whose on- 
shell matrix elements are the transition amplitudes, are 
obtained by solving the three-body AGS integral equa- 
tion 0; 01 that reads 

C/'^" = SpaGo' + E U^^Got^S-^a , (5) 
7 

or 

U"" = SpaGo^ + E Sp^t^GoU^" , (6) 

7 

where Spa = 1 — <5/3a- The pair transition operator is 

tj = + v~^Gotj , (7) 
where Go is the free resolvent 

Go = {E + lO - Ho)-\ (8) 

and E the total energy of the three-particle system in the 
center of mass (cm.) frame. For breakup {(3 — in the 
final state) one has 

C/0"^G^i+^i^GoC/'^", (9) 

7 

where U^" is obtained from the solution of Eq. ^ with 
a, j3,j = (1, 2, 3). The scattering amplitudes are the ma- 
trix elements of U^" calculated between initial and final 
states that are eigenstates of the corresponding channel 
Hamiltonian Ha {Hp) with the same energy eigenvalue E. 
For breakup the final state is the product of two plane 
waves corresponding to the relative motion of three free 
particles that may be expressed in any of the relative Ja- 
cobi variables. In the latter case the contribution of the 
Gn ^ term is zero. 
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FIG. 1: Single scattering diagrams for breakup in the Fad- 
deev/AGS scattering framework. 



The solution of the Faddeev/AGS equations can be 
found by iteration leading to 

jjt)a _ t^S^a + E] ^1 5^ Go5^^t(5^a 

7 7 C 




+ •■• , (10) 

where the series is summed up by the Fade method [iq| . 
The successive terms of this series can be considered as 
first order (single scattering), second order (double scat- 
tering) and so on in the transition operators. For the 
breakup process where particle 1 scatters from pair (23) 
into the breakup channel (1,2,3), the breakup series in 
single scattering is represented diagrammatically in Fig.[T] 
where the upper particle is taken as particle 1. 

In our calculations Eqs. l(6l fT0l) are solved exactly in 
momentum space after partial wave decomposition and 
discretization of all momentum variables. We include 
the nuclear interaction between all three pairs, and 
the Coulomb interaction between the proton and ^'^Be, 
following the technical developments implemented in 
Refs. [nl, for proton-deuteron and a-deuteron elastic 
scattering and breakup that were also used in Refs. [1, Q 
to study p-^^Be elastic scattering and breakup. 



III. THE PHYSICAL CONTENT OF THE 
SINGLE SCATTERING TERM AT NP QFS 

Let us consider a reaction where the proton (particle 
1) scatters from a bound pair of particles (23) into the 
breakup channel (1,2,3); particle 2 is a valence neutron 
and particle 3 the ^'^Be core. Actual experiments involv- 
ing halo nuclei are performed in inverse kinematics with 
the scattering of a radioactive halo beam from a stable 
target. The results discussed here are however indepen- 
dent of the kinematics that is used. 

We have shown Q that at sufficiently high energies 
and for some suitable kinematics configurations the sin- 
gle scattering term provides a reasonable approximation 
to the multiple scattering series. In addition, the single 
scattering component where the projectile strikes the va- 
lence particle (first diagram in Fig. [l]) is dominant in the 
case of p("Be,i*'Be n)p 0. Therefore, we analyze here 
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in detail the physical content of the two components of 
the single scattering term involving both the scattering 
from the struck neutron and the ^^Be core. 

In this section we shall use capital letters for momenta 
in the LAB frame and lower case letters for the cm. 
frame. 

In the case of the single scattering from the struck 
neutron, the transition amplitude from an initial state 
= |ki(/)23) to a final state jk'j^kjkg) is given by the 
matrix elements of the transition operator = ti2{uji2) 
with relative pair (1,2) energy ^12. As usual ki is the ini- 
tial relative momentum between particle 1 and the cm. 
of pair (23). In the final state, k2,k2, and kg are the fi- 
nal momenta for the three free particles in the three-body 
cm. frame where k^ -I- k2 -I- kg = 0. 

The breakup transition amplitude, corresponding to 
the single scattering approximation (SSA), may be writ- 
ten as 

{KKK\t3\'4'l) ^ (qi2l^l2(wi2)|qi2)023(q23) , (11) 



where 



q'i2 - K 



qi2 = ki 



mi 



mi 



K , (12) 



M 



q23 



UJ12 = E 



-k' 



12 



m3 



k^, , (13) 



M- 



ki 
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(14) 



9l2 



2/^3(12) 2/X12 



(15) 



with Mij — rrii 



+ 1Tlj, flij = Ti~~^ /^i(jfc) 



M 



and 

M = rrii + rrij + m^. Therefore the two-body t-matrix 
ti2 is calculated at half the energy shell. 

At the np quasifree scattering (QFS) kinematical con- 
dition particle 3 (the ^'^Be core) acts as a spectator which 
means that in the Lab frame K3 = corresponding to 
kg = —j^'ki in the cm. frame. Therefore, according to 
Eq. (O 



[q23] 



QFS _ 



and 



W12 



I QFS 



M 1712 

M23 Mr. 



IQFS 



2^1 



(16) 



(17) 



in the limit of zero binding energy for the pair 23 as 
demonstrated in the Appendix A. Only in such zero 
binding limit and in the QFS kinematical condition does 
the t-matrix ^12(^12) gets calculated on the energy shell. 
Away from the QFS kinematical condition the ii2(wi2) 
matrix elements become half-on-shell and the total tran- 
sition amplitude probes the non zero relative momentum 



components of the bound pair. Higher order multiple 
scattering terms necessarily probe off-energy-shell effects 
and higher relative momentum components, even at the 
QFS conditions. 

We now consider the SSA component where the proton 
scatters from the core. The transition amplitude is then 
given by 

(k'ik^mt2|V'i) = (q3l|il3(wi3)|q3l)'^23(q23), (18) 



with q(,i = -^kf. 



q23 = k^ -f 



k' 



2 ^ '^i: q3i 

and Wi3 = E 



' Mi3 *^2 



ki, 



M23 J- > - - ^^2(13) 

In the np QFS limit where kg = — kj - 
we get 



k; = 



k2 = ^k'l + TT^ki 



23 



"^3 1,, 



(19) 



leading to 



[q23]^^' = ki - k'l = k 



where k is the momentum transfer of particle 1, the pro- 
jectile. Therefore, in the case of the projectile scattering 
from the core, the wave function of the target nucleus is 
probed at nonzero relative momentum. For the case of 
the reaction p(^^Be,^''Be n)p studied in the valence 
neutron is bound to ^°Be in S-wave whose wave func- 
tion momentum distribution is sharply peaked around 
zero momentum. Thus, in this case, the wave function, 
when probed at larger momentum ^23, is already very 
small, leading to a scattering contribution from the core 
that is relatively small. For other orbital angular mo- 
mentum states such as P- or D-waves, the wave function 
becomes nonnegligible at larger (723 and the contribution 
from projectile-core scattering may become significant. 
Therefore standard DWIA calculations which only take 
into account the contribution from the scattering of the 
struck particle may become inadequate in the case of 
struck particles bound with nonzero orbital angular mo- 
mentum. We shall return to this point later. 



IV. THE BREAKUP OBSERVABLES IN 
INVERSE KINEMATICS 

We now consider the breakup of a radioactive beam 
involving a two-body halo nucleus assumed to be well 
described by a core and a valence neutron; the halo nu- 
cleus collides with a proton target leading to three free 
particles in the final state. This final state is described 
in terms of 9 kinematical variables. Momentum and en- 
ergy conservation reduces this number to 5 independent 
variables. With the recent developments at the radioac- 
tive beam facilities it is now possible to measure fivefold 
fully exclusive observables. By further integration upon 
the variables of the emitted particles semi-inclusive as 
well as inclusive breakup observables can be measured 
although with loss of physical information. 
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In actual experiments performed in inverse kinematics 
it is the halo core which is measured either directly or 
by reconstruction from the other emitted fragments. We 
therefore chose the Jacobi momenta 

m„K„ - m„K„ 

P = \ , 

nip + nin 

(mp + ?7i„)Kc - rnc(Kp + K„) 

q = — , (20) 

with Kp, K„, Kc {nip, m„, mc) being the LAB momenta 
(masses) of the proton, valence neutron and core in the 
exit system, and M = rup + rUn + mc- The breakup dif- 
ferential cross section is calculated from the on-shell ma- 
trix elements of the AGS operators, — (qplC/*^" |?/'a) 
where particle a is the spectator in the initial state (in 
our case a is the proton) and the Jacobi momenta in the 
final state satisfy the on-shell relations 



with the reduced masses 



? 2 

2/i 2/1 







(21) 



nin + mp 
mc{m„ + TUp) 
M 



(22) 



Fully exclusive observables 

The fully exclusive breakup observables are measured 
in the LAB system. The kinematic configuration of three- 
body breakup is characterized by the polar and azimuthal 
angles fli = {6i,(j)i) of the two detected particles as in 
Fig. [2] which we assume to be the core C and the valence 
neutron n. 

In order to calculate this observable we begin by writ- 
ing the general expression in terms of the on-shell matrix 
elements of the AGS operators an including momentum 
and energy conservation 



d^cr 
dKdKc 



{2tt 



4 m„ -I- mc 



I dKpKldK^\T''^\^S{Ki^AB - K„ - Kp - Kc)S{Ei^ab + e 

(Klab - K„ - Kc)2 



R'2 



2mp 2m,j, 2mc ' 



^^^^, m +mc I K2dK^^j,0a^HiE^AB + e 



= (27r) — — mpm„ > 



2mr. 



2to„ 2mc 



IT' 



Oct 1 2 



(to„ + mp)Kn ~ to„(Klab - Kc) • K„| 



(23) 



where the sum on i involves the momenta Kn given by 
the zero's of the argument of the energy conserving 5- 
function 



E, 



LAB 



(Klab - K„ - Vicf K 
2mr, 2mn 



2mc 



, 



(24) 



Klab is the beam momentum in the Lab frame and E'lab 
the corresponding energy. To arrive to Eq. (|23| we used 
the property of the i5-function. 



%w) = E 



6{x - Xj) 



(25) 



often replaces the dependence on the energy Ec by the 
arclength S related to the LAB energies En and Ec of 
the two detected particles as S" = J„^ ^5* with 



dS 



(26) 



By differentiating the argument of the energy conserving 
(5-function with respect to Kc and taking Kn — Kn{Kc) 
we obtain 



where the the zeros of the function g{x). Eq. ([24]) 

defines an ellipse in the Kn — Kc plane. The points 
lying on that ellipse or on the corresponding curve in 
the En — Ec plane comprise the kinematically allowed 
S-curve on which the physically accessible events have 
to lie 0. The phase space factor in Eq. l(23|) has the 
disadvantage that it diverges at the Kc values at which 
the ellipse has an infinite derivative. As a result one 



dKn _ {mp -I- mc)Kc - toc(Klab - K„) ■ Kc 
dKc {mp + mn)Kn - m„(KLAB - Kc) • K„ mc 



(27) 



Therefore one can measure the fully exclusive fivefold 
differential breakup cross section d^a/dflndflcdS where, 
from Eq. l(23|) . we get 
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dilndflcdS 



{2tt) — — mpirinmcKc 

^LAB 



[ I (to„ + mp)ii:„ - to„(Klab - Kc) • K„ 



\ TUnKc dKc J 



(28) 




FIG. 2: (Color online) Kinematic angles for breakup in inverse 
kinematics in the LAB frame 



with dKn/dKc given by Eq. (|27|) . The sum on i disap- 
pears in Eq. (|28l) relative to Eq. (|23|) . because there is a 
one to one correspondence between £"„ and S which may 
not exist between En and Ec- 



Semi-inclusive observables 

The semi-exclusive differential cross section 
d^a/d^lcdEc may be obtained from the fivefold 
differential cross section ((23|) by integrating over the 
angles of the emitted valence neutron. However, it is 
more convenient to start from the fivefold differential 
cross section in the cm. frame 



d^p d^q 



-mc 



Klab 



Using the property l|25p of the i5-function one gets 



d'^p d^q 



(27r)4!^^^±!^|T0"|Vj2^i?-^g2 . (30) 



From Eq. ^ , using /Ctot — /Cn + ^p 
bering that in the chosen Jacobi set q 
one gets 



-Kc and remem- 



-TOT 



= (2-)'^^^^|^-c^c / d'p\T"'\'^^^2^,E -^{)Ch + §^loT - • /Ctot) , (31) 



where the calligraphic momenta may be LAB or cm. 
momenta depending on /Ctot = -fs^LAB or /Ctot = 0. 
Clearly in the integrand of this equation the square 
root must be positive definite and, therefore, when the 
LAB energy of the emitted core exceeds a certain limit, 
E^^^ (max) , which depends on the angle of the emitted 
core, then the semi-exclusive cross section should van- 
ish. In order to find the rate at which the cross section 
vanishes one should take the derivative of the differential 
cross section with respect to the energy of the core. It 
follows that at _E^'^^(max) this derivative is infinite and 
therefore the semi-inclusive cross section should exhibit 
almost a sharp cutoff at the maximum value of the core 
laboratory energy unless the breakup amplitude is very 
small close to the maximum energy. Hence the energy 
behavior of the semi-inclusive cross section at high ener- 
gies follows from a delicate interplay between the phase 



space and the scattering amplitude which should thus be 
accurately calculated. 



Inclusive observables 

As shown in the Appendix B, from this semi- inclusive 
cross section one can calculate the inclusive perpendicu- 
lar momentum distribution for one of the detected parti- 
cles, say the core, as 



dRP / ^3j^ c 



(32) 



J -co ^ J^^C 

and the inclusive transverse momentum distribution as 



da 



= 2 



+ 00 /' + OC 



^dK^cdKl (33) 
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with 

_ 1 (Pa 
(PKc ~ mcKc dVLcdEc ' 

V. THE PAIR INTERACTIONS 

In this work we address the scattering of the (^'^Be +n) 
system by a proton target. Before solving Faddeev/AGS 
equations we need to define each pair interaction, that is, 
the p-n, p-^^Be and n-^^Be interactions. 

For the p-n we take the reaUstic nucleon-nucleon CD 
Bonn potential 

For the potential between the proton and ^^Be core 
we use a phenomenological optical model with param- 
eters taken from the Watson global optical potential 
parametrization [1, [lB| • The energy dependent parame- 
ters of the optical potential are taken at the proton labo- 
ratory energy of the p-^^'Be reaction in direct kinematics. 

The interaction between the valence neutron and the 
^^Be core in ^"^Be depends on the orbital angular config- 
uration of the valence nucleon for the ^^Be + n cluster 
system. Since the total angular momentum of the core 
and valence neutron are coupled to ^'*Be(0"'"), it means 
that the the total angular momentum of the valence nu- 
cleon is identical to the ^'^Be core. Little is known about 
this system with exception of a possible 5/2+ (D-wave) 
resonance with a relative energy i?roi(^^Be+n)^2 MeV 
above the neutron breakup threshold. 

In this work we take the spectroscopic study of ^^Be 
obtained from the proton-induced reaction on ^''Be at 
69 MeV/u in inverse kinematics performed by the invari- 
ant mass method at RIKEN The ^^Be resonances of 
1/2- (P-wave) with E,c\ = 0.45 MeV, 1/2+ (S-wave) with 
£:,.ci = 1.17 MeV, and 5/2+ (D-wave) with E.^i = 2.34 
MeV were assigned through the analysis of the transverse 
momentum distributions of the outgoing ^'^Be system. 
Following this work we take three possible single-particle 
configurations for ^''Be(0+) 

P^Be) = li3Be(l/2-) n(lpl/2)), (35) 
I^^Be) = |i3Be(l/2+) ® n(2sl/2)), (36) 
I^^Be) = |"Be(5/2+) ® n(ld5/2)), (37) 

with binding energy e — E^ci + S2n, S2n — 1-26 MeV 
being the two neutron separation energy of ^^Be. 

The interaction between the valence neutron and the 
^^Be core is taken to be of the form 

Vir) = -VJir,Ro,ao) (38) 

where /(r, R, a) is the usual Woods-Saxon form factor 

fir, R, a) = 1/{1 + cxp[(r - i?)/a]}, (39) 

and Ri = ViA^ . The potential parameters for each single- 
particle configuration are listed in Table H] and the corre- 
sponding wave functions shown in Fig. [3l 



TABLE I: Parameters of the n- Be interaction plus ro = 1.2 
fm and ao = 0.6 fm. 



State 


K(MeV) 


lpl/2 
2sl/2 
ld5/2 


32.922 
61.667 
68.980 
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-51 ' ^ ' 1 ' 1 ' 1 

12 3 4 

P (fm"') 

FIG. 3: (Color online) Relative wave function of the valence 
neutron and ^^Be core in momentum space for S-(dark solid 
line), P-(dark dashed line) and D-wave bound states (light 
solid line) 

In order to access the contribution of other partial 
waves in the n-^^Be interaction beyond the one that is 
responsible for the single particle configuration leading 
to the 0+ ground state of ^"^Be, we take in those partial 
waves the potential corresponding to 14 = 32.922 MeV 
in Table [H Therefore 0+ partial waves are driven by the 
potentials in Table [H depending on the choice of ^^Be 
resonance we consider; all other partial waves are driven 
by the weakest potential of all three for no other reason 
than lack of a more enlightened choice. 

VI. RESULTS 

In the solution of the Faddeev/AGS equations we in- 
clude n-p partial waves with relative orbital angular mo- 
mentum inp < 8, and p-^^Be with L < 10. For the n-^^Be 
we start by including just the 0+ partial wave using the 
interactions in Table H] depending on the specific choice 
of spin and parity we consider for ^'^Be, that is, l/2~, 
1/2+ , or 5/2+. Three-body total angular momentum 
is included up to 100 and in all calculations unphysical 
bound states are properly removed j§] keeping a single 
bound state of ^^Be in 0+ partial wave. 

We present in this section the calculated semi-inclusive 
breakup observables around QFS and the inclusive mo- 
mentum distributions of the ^^Be system. 

In inverse kinematics the exact np QFS kinematic 
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conditions correspond to a ^'^Be LAB momentum of 
Kc = j^Kc„, that is, taking ^''Be in the ZZ direc- 
tion we ha^e (iS^AB^ ^lab^ ^lab^qfs ^ (^^^^ o,0), 
where the index C denotes ^^Be, Cn denotes ^''Be and 
Mcn = fnc + "^n- For ^^Be scattering from a proton at 
69 MeV/u this corresponds to (897 MeV,0,0). 

In Fig. [4] we show the semi-inclusive cross section for 
the breakup pC'^Be,^^Be)np at 69 MeV/u using the SSA. 
We separately consider the scattering from the valence 
neutron, the scattering from the core, and the scattering 
from both the valence and the core. As shown in Sec. Ill, 
the cross section calculated using SSA from the struck 
particle probes the bound state wave function at relative 
momentum (723 = in the np QFS kinematic condition. 
As shown in the Fig. [3] this corresponds to the point 
where the relative wave function is sharply peaked in the 
case of S-wave, but zero in the case of P- or D- waves. 
At Ec = 897 MeV, as one moves away from the QFS 
point, that is, when 0c increases from zero, the cross 
section starts to probe the wave function at larger relative 
momentum. For the S-wave case where the wave function 
is very narrow in momentum space, this means that the 
cross section will decrease very rapidly as one increases 
the core scattering angle, as seen clearly in Fig.[4l For the 
case of a P-wave, it starts to probe nonvanishing values 
of the wave function reaching steeply its peak. Likewise 
for the D-wave as one moves away from the QFS point it 
starts to probe nonvanishing values of the wave function 
reaching less steeply its peak as also shown in Fig. [H 

On the other hand the cross section calculated using 
SSA from the core probes the bound state wave function 
at relative momentum different from zero. This means 
that the contribution from the scattering from the core is 
small for the S-wave. As shown in the graph at this en- 
ergy the scattering from the core gives an important con- 
tribution to the single scattering term in the case where 
the struck particle is bound to the core in P- or D- waves. 

We would like to point out that the SSA from the 
core is neglected in standard DWIA approaches which 
are therefore inadequate for nonzero relative angular mo- 
mentum wave functions of the bound pair. 

In Fig. [5] we show the semi-inclusive cross section for 
the breakup p(^^Be,^^Be)np at 69 MeV/u using full mul- 
tiple scattering Faddeev/AGS calculations. The most 
striking feature from the graphs is the relative order of 
magnitude of the calculated observables for each single- 
particle configuration which is carried over from the SSA 
calculations shown in Fig. IH The breakup observable 
calculated with the S-wave single particle state is one or- 
der (two-orders) of magnitude larger that the breakup 
observable calculated with the P-wave (D-wave). There- 
fore the magnitude of the semi-inclusive cross section in 
quasifree conditions is a clear signature of the angular 
momentum of the valence nucleon. Comparing Figs. [4] 
and [5] we conclude that the SSA clearly overestimates 
the full multiple scattering results in all configurations 
at this energy regime. Therefore care should be taken 
into account when using truncated multiple scattering 



frameworks such as DWIA as pointed out in Ref. 

In order to assess the effect of introducing higher par- 
tial waves in the n-^'^Be interaction we follow the pre- 
scription mentioned at the end of Section |V] and use the 
potential corresponding to Vc = 32.922 MeV in Table H] 
in all partial waves other than 0"''. For that matter we 
include n-^'^Be relative angular momenta ^ < 3. The re- 
sults for the semi-inclusive cross section obtained from 
the full multiple scattering Faddeev/AGS calculation are 
shown in Fig. [6] and do not differ significantly from the 
results in Fig. [H This indicates that the effect of intro- 
ducing the coupling to the other n-^^Be partial waves is 
small, and does not change the conclusions. 

In Fig. [7] we show the inclusive transverse momentum 
distributions of the ^'^Be = (^^Be+n) system in the l/2~, 
1/2+ , and 5/2+ states at 69 MeV/u obtained from the full 
multiple scattering Faddeev/AGS calculations. These re- 
sults show that the shape of the inclusive transverse mo- 
mentum distribution does not provide a clear signature 
for the spin of the ^^Be resonance, but instead, its mag- 
nitude at the peak does, as long as one has full control 
of the three-body dynamics. Nevertheless the magnitude 
of the semi-inclusive cross section is far more sensitive to 
the orbital angular momentum of the (^^Be+n) 0+ bound 
state than the inclusive transverse momentum distribu- 
tion of ^^Be. 



VII. CONCLUSIONS 

We have performed full Faddeev-type calculations for 
one-neutron knockout reaction of ^"^Be on proton target 
at 69 MeV/u incident energy. These results were com- 
pared with those corresponding to the single scattering 
approximation. Inclusive transverse momentum distribu- 
tion observables for the outgoing (^^Be+n) system were 
calculated. Semi-inclusive differential cross sections at 
the np QFS kinematical conditions were also presented. 

In this work we have considered the 0+ ground state 
of ^"^Be as a single particle state made up of (^'^Be+n) 
bound in P-, S-, or D-wave depending on the ^'^Be spin 
being 1/2^, 1/2+ , or 5/2+, respectively. 

We have found that the single scattering contribution 
from the core is very significant in the case of P- or D- 
waves. In any case the single scattering is a bad approx- 
imation of the full results at this energy. Thus, DWIA 
approaches which take into account the contribution be- 
tween the proton and the struck valence neutron, are 
inadequate in this case. Higher order multiple scattering 
contributions need to be taken into account as done in the 
full Faddeev/AGS approach. The semi-inclusive breakup 
cross section resulting from the full Faddeev/AGS calcu- 
lation with a (^•^Be+n) S-wave bound state is one order 
(two-orders) of magnitude larger than in the case of the 
P-wave (D-wave) bound state. Therefore the magnitude 
of the semi-inclusive cross section in QFS conditions is a 
clear signature of the angular momentum of the valence 
nucleon if we can have a good control of the underlying 
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FIG. 4: (Color online) Semi-inclusive cross section for the breakup Be,^ Be)np at 69 MeV/u using the SSA contribution 
for 6c ~ 0° (dark thick solid line), 0c = 0.4° (dashed line), dc = 0.8° (dark thin solid line), 9c = 1-2° (dashed dotted line), 
and 9c = 1.4° (light solid line). 



three-body dynamics. 



leading to 



APPENDIX A: 

In SSA, the initial relative momentum betvireen parti- 
cles 1 and 2 is 



qi2 



m2 



ki 



mi 



(-ki-k^). 



(A.l) 



In np QFS kinematical conditions kg = — ki „^J!|?,„^ 
get 



[qi2] 



QFS 



Mm2 
M12M23 



(A.2) 



The two-body energy W12, in the limit of zero binding 
energy for the pair (23), is given by 



U12 = E - 



2M3( 



12) 



2^1(23) 2^3(12) 

Under the QFS condition 

QFS _ ^1 ^1 



(A.3) 



W12 



E 



2^1(23) 2/^3(12) M23 



2miM23^Mi2 
Mm2 



M23M: 



12 



(A.4) 



W12 



QFS 



2a^(12) 



(A.5) 



Thus, in the single scattering term the matrix elements 
of the transition operator ti2(wi2) are on the energy shell 
in the limit of zero binding for pair (23). 



APPENDIX B: 

Let us consider the breakup reaction p((nC),C)np. In 
this section we give the formulae for the momentum dis- 
tributions of the detected core C. The kinematics of 
particle C can be defined in terms of its momentum 
{Kx, Ky, Kz) or, alternatively, by its LAB energy and an- 
gular variables {E, 51) where, in this section, we drop the 
index of the core for simplification. The semi-inclusive 
cross section is given by 



1 (Pa 



<PK mK dEdVL 



(B.l) 



with m the mass of the core. The Cartesian compo- 
nents of the momentum of the particle can be expressed 
in terms of its spherical components {Kr^,Ky,Kz) = 
(K sinO cos(f>, if sin ^ sin 0, Kcos9). Alternatively one 
may define the set of cylindrical momentum coordinates 



KlMn-\KylKx),K,). One 
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FIG. 7: (Color online) Transverse momentum distributions of 
the ^^Be = (^^Be+n) system in the p("Be,^^Be)np reaction 
at 69 MeV/u. Results of the full Faddeev/AGS calculations 
with S-, P- and D-wave (^^Be+n) bound states are shown. 



FIG. 5: (Color online) Semi-inclusive cross section for the can write then 
breakup p(^''Be,^^Be)np at 69 MeV/u using the full multiple 
scattering calculations for Qc ~ 0° (dark thick solid line), 
dc = 0.4° (dashed line), Qc = 0.8° (dark thin solid line), 
61c = 1.2° (dashed dotted line), and Qc = 1.4° (light solid 
line). The n-^'^Be interaction is restricted to the 0^ partial 
wave. 



d^K = dK^dKydK^ = d^KpdK^ = KpdKpd(f>dKJyB.2) 



Full Faddeev-all partial waves 



4x10 



2x10 



> 



4x10 
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a 

-n"^ SXIO' 



4x10 
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Thus 



d^Ko 



+ 00 



dK, 



(B.3) 



From this double cross section we can calculate the in- 
clusive perpendicular momentum distribution 



= 2t:K„ / ^^dK^ 



dK, 



d^K 



(B.4) 



and the inclusive transverse momentum distribution 



(B.5) 



FIG. 6: (Color online) Semi-inclusive cross section for the 
breakup p(^'*Be,^^Be)np at 69 MeV/u using the full multiple 
scattering calculations for 6c = 0° (dark thick solid line), 
dc = 0.4° (dashed line), Oc = 0.8° (dark thin solid line), 
Oc = 1.2° (dashed dotted line), and Oc = 1.4° (light solid 
line). The n-^'^Be interaction is included in all partial waves 
(see text). 
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